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ABSTRACT 

Gap formation in a gas disk triggered by disk-planet tidal interaction is considered. 
Density waves launched by the planet are assumed to be damped as a result of their 
nonlinear evolution leading to shock formation and its subsequent dissipation. As a 
consequence wave angular momentum is transferred to the disk, leading to evolution 
of its surface density. Planetary migration is an important ingredient of the theory; 
effects of the planet-induced surface density perturbations on the migration speed are 
considered. A gap is assumed to form when a stationary solution for the surface density 
profile is no longer possible in the frame of reference migrating with the planet. An 
analytical limit on the planetary mass necessary to open a gap in an inviscid disk is 
derived. The critical mass turns out to be smaller than mass Mi for which planetary 
Hill's radius equals disk scaleheight by a factor of at least Q^^"^ {Q is the Toomre stability 
parameter) depending on the strength of the migration feedback. In viscous disks the 
critical planetary mass could vary from ~ O.lMi to Mi, depending on the disk viscosity. 
This implies that a gap could be formed by a planet with mass (1 — 10) depending 
on the disk aspect ratio, viscosity, and planet's location in the nebula. 

Subject headings: planets and satellites: general — solar system: formation — (stars:) 
planetary systems 



1. Introduction. 



The recent discovery of extrasolar planets on orbits very close to their parent stars (Mayor & 
Queloz 1995; Marcy et al. 2000; Vogt et al. 2000; Butler 2001) has raised a number of questions 
about the formation mechanisms of such systems. These close-in planetary companions are all 
presumed to be gas giants, typically with masses of the order of Jupiter mass Mj 2 y. lO^'^ g. 
It is very unlikely that such planets were formed at their present locations (Boss 1995): current 
theories (Mizuno 1980; Bodenheimer & Pollack 1986) predict that giant planets were formed by 
gas accretion onto massive (~ 15 M^) rocky core which themselves are the result of accumulation 
of a large number of icy planetesimals. The most favorable conditions for this process exist beyond 
the so-called "snow line" (Hayashi 1981; Sasselov & Lecar 2000) which is estimated to lie several 
AU from the star, far larger than the actual orbital radii of the known extrasolar planets. 
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The most popular theory explaining this paradox presumes that giant planets were indeed 
formed far outside their present locations but then migrated inwards as a result of tidal interaction 
with a gaseous nebula (Ward 1997a). If a planet is fully immersed in a gas disk then the migration 
process works as follows: planetary gravitational torques produce density waves carrying angular 
momentum in the surrounding gas disk. For the planet to migrate inwards (outwards) two condi- 
tions must be fulfilled: (1) interaction with the outer (inner) part of the disk should be stronger 
than with the inner (outer) one, and (2), the waves must not return to the planet. In a finite disk 
the only way to fulfill the second condition is for the waves to dissipate so that their angular mo- 
mentum gets transferred to the disk flow. It was demonstrated by Ward (1986) that in Keplerian 
disks conditions are usually such that the planet tends to migrate inwards. The typical timescale 
of this process, ~ 10^ yr for a Jupiter-mass planet, is very short compared with the lifetime of the 
nebula itself (10^ — 10^ yr). This short timescale presents a significant problem for the migration 
scenario because the migrating planet is likely to drift right into its parent star in very short time. 

Various mechanisms to stop migration have been proposed: formation of a cavity in the inner 
disk by magnetospheric activity of the central star (Shu et al. 1994), resonant interaction with 
another planet in the system (Masset & Sncllgrove 2001), spin-orbital interaction with the central 
star (Trilling et al. 1998), etc. On the other hand, there is a process almost inevitable in the course 
of giant planet formation — gap formation in the disk near the planet — which can effectively slow 
down planet migration. When the planet is not very massive its gravity cannot strongly affect the 
disk surface density distribution in its vicinity and strong interaction with the gas disk through 
Lindblad resonances (Goldreich & Tremaine 1980, hereafter GT80) leads to the rapid migration 
described above (so called "type I migration", Ward 1997a). As the planetary mass grows the 
drift speed increases but at some point the strength of the torque exerted on the surrounding gas 
becomes so large that planet simply repels the gas and a gap forms. This strongly diminishes the 
tidal interaction because it is usually dominated by high-order Lindblad resonances which now lie 
inside the gap. Thus the orbital evolution of the planet becomes tied to the evolution of the disk 
and it migrates on the viscous timescale of the disk which could be much longer (10^ — 10^ yr or 
even longer, depending on the viscosity in the disk). This stage is called "type II migration" and 
its existence might help to explain the survival of planetary systems in the course of their orbital 
evolution or at least significantly alleviate this problem (Ward 1997). Gap formation also provides 
a reason for the existence of a maximum mass of extrasolar giant planets (Nelson et al. 2000). 

In view of all this a very important question arises: how massive should a planet be in order 
to open a gap? The answer depends not only on the conditions in the nebula (surface density, 
temperature, viscosity) but also on the dissipation mechanism of planet-induced density waves and 
the mobility of the planet. If the viscosity in the disk is absent and the planet is not drifting there 
is no mechanism which could oppose tidal torques and a gap is opened by an arbitrarily small 
perturbcr. This conclusion changes if planetary migration is taken into account sclf-consistcntly 
because a low-mass planet is able to drift through the gap before gap fully forms (Ward &: Hourigan 
1989, hereafter WH89). The planet must slow down its migration somehow for the gap to be opened. 
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In their study WH89 demonstrated that the disk surface density is enhanced in front of the 
drifting planet and is reduced behind it. This effect diminishes the difference between the torques 
produced by the planet in the disk outside and inside its orbit and acts to slow down the migration. 
For some mass of the planet its steady drift becomes impossible, the planet halts, and gap formation 
ensues. Thus, in this picture there is a minimum mass for gap opening even in an inviscid disk. 
Assuming that wave dissipation is a very rapid process and that angular momentum is immediately 
transferred to the disk WH89 estimated this critical mass to be 

h 

M2 ~ -Mf, where M/ = S/i^. (1) 

Here S is a surface density of the disk, h = c/fl is its geometric thickness (here c is the sound 
speed and Q is the angular orbital frequency), and r is the distance from the central star. Typical 
protosolar nebular parameters were summarized by Hayashi (1980) in the form of the minimum 
mass Solar nebula (MMSN) model: 

q / r \-3/2 1 / r \-V4 

So(r) = 1700 g cm-3 , co(r) = 1.2 km g-^ (^) . (2) 

Using (2) we obtain that M/ « 6 x lO^^g « 0.1 and M2 ~ 3 x lO^^ g at 1 AU which is about 
the mass of Mercury. Of course, if the disk possesses nonzero viscosity this mass would increase 
because viscous diffusion opposes gap formation. 

Lin &; Papaloizou (1993, hereafter LP93) present a different point of view on the planetary 
mass required to open a gap. They also assume damping of the density waves to be local but require 
mass of the planet to be high enough for the waves to be strongly nonlinear from the beginning and 
shock immediately transferring angular momentum to the fluid. They showed that this minimum 
mass corresponds to the case when the Hill's radius of the planet Rh = r{Mp/M^,y^'^ (Mp is the 
mass of the planet, is the mass of the central star) is comparable to the vertical disk scaleheight 
yielding 

2c^ / r \3/4 

= SfiO ^ " (lAu) 
as the relevant mass for opening a gap [the numerical estimate in equation (3) is made for MMSN 
parameters given in equation (2)]. Another way to look at this restriction is to notice that at 
Mp ~ Ml the pressure gradient in the disk in the vicinity of the planet becomes so high that 
epicyclic frequency becomes imaginary and Rayleigh's instability develops. Thus, this restriction 
could be considered as an upper limit on the mass of the planet that does not open a gap. 

Values of the gap opening mass estimated by these two approaches differ by a huge factor. 
Indeed: 

m ^ 2nj^ ^ 27r 

Mf 3 ttGS 3 ^' ^ ^ 

where Q is the Toomre stability parameter (Binney & Tremaine 1987) and k is the epicyclic 
frequency {k = for Keplerian rotation law); using MMSN parameters given by equation (2) we 
obtain that Q ^ 70 and h/r 0.04 at 1 AU meaning that M2 ~ {h/r)Q-'^Mi 5 x IQ-^Mi. 
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Ward & Hourigan (WH89) assume that the damping mechanism is independent of the plan- 
etary mass and essentially linear, such as viscous dissipation (Takcuchi ct al. 1996) or radiative 
damping (Cassen &: Woolum 1996). However these mechanisms arc probably ineffective in cold, 
weakly ionized and optically thick systems such as protoplanetary disks (Hawley & Stone 1998; 
Goodman & Rafikov 2001, hereafter GROl); in this case tides raised by low-mass planets could 
propagate much further than just a fraction of disk scale length and gap opening requires signifi- 
cantly more massive perturbers (Ward 1997a). 

Requiring waves to be strongly nonlinear and damp immediately as a necessary condition for 
gap formation is also probably too radical. Indeed, density waves produced by the planet can still 
evolve due to weak nonlinearity and are able to form a weak shock even if the planet is less massive 
than Ml (GROl). This mechanism can lead to a gap formation by lower mass planets than LP93 
assumed. Clearly to obtain a reliable estimate of the critical planetary mass one must use a realistic 
wave damping prescription. 

Goodman & Rafikov (GROl) have considered nonlinear evolution of the density waves produced 
by low-mass planets in two-dimensional disks using the shearing sheet approximation and assuming 
the background surface density and sound speed to be constant. They have found that a shock is 
formed quite rapidly (depending on Mp/Mi, typically after travelling several disk scaleheights from 
the planet but still not immediately) because the radial wavelength of the perturbation constantly 
decreases while its amplitude increases (as a consequence of angular momentum flux conservation) 
thus facilitating wavefront breaking. After wave shocks its angular momentum is gradually trans- 
ferred to the disk fluid leading to surface density evolution (thus violating the constant surface 
density assumption). Rafikov (2001) generalized this analysis by taking into consideration the ef- 
fects of spatial variations of the surface density and sound speed in the disk as well as the cylindrical 
geometry of the problem. Realistic prescription of the global angular momentum dissipation was 
provided under the condition that the surface density and sound speed vary on scales larger than 
the wavelength of the perturbation. 

In this paper we study the criterion for gap opening by planetary tides by determining un- 
der which conditions a steady-state solution for the disk surface density perturbation is no longer 
possible (method used in WH89). Planet migration is taken into account self-consistently. Tidal 
perturbations are supposed to be damped nonlocally by weak nonlinearity and their angular mo- 
mentum is assumed to be transferred to the disk, as described quantitatively in Rafikov (2001). For 
the gap opening the assumption of varying surface density in Rafikov (2001) is the most important 
one because it allows us to solve the problem self-consistently (angular momentum transfer depends 
on how the surface density is distributed radially). 

After developing the general analytical apparatus in §2 wc first explore the case of an inviscid 
disk in §3.1 to highlight the most important physical mechanisms relevant for gap formation. Then 
we generalize consideration to include the disk viscosity in §3.2. Finally, in §4 we discuss our results. 
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2. Basic equations 



We study the surface density evolution of a gas disk driven by planetary gravitational per- 
turbations. In our model a planet of mass Mp moves in a disk on a circular orbit with radius r^. 
The background or initial surface density and sound speed 'So{r) and co(r) are assumed to vary 
independently on scales ~ rp. The values of the surface density, sound speed, orbital frequency, 
and disk vertical scaleheight at planet's location are T,p,Cp,Qp, and hp respectively. We have sum- 
marized brief descriptions of and references to definitions of the most important variables we use 
in Table 1. 

We are only interested in the radial evolution of the surface density and neglect azimuthal 
variations. Then the general time-dependent equations of disk surface density evolution are the 
continuity and angular momentum equations (Pringle 1981): 



r— + —{r^Vr) = 0, 



dt dr 



dt 



dr 



J_ fdF_ dG 
2tt V dr dr 



G = -27rr'VS 



dr ' 



(5) 

(6) 



where S is the disk surface density, Vr is a fluid radial velocity, F is the (azimuthally averaged) 
angular momentum flux transferred to the disk fluid from the planetary tidal perturbations, and 
G is the usual viscous angular momentum flux [v is a kinematic viscosity). 



Prom this system one can easily find that 

1 



27rrS 



n -1 



(7) 



Substituting this result back into equation (5) one obtains a self-consistent equation for the surface 
density evolution: 



as 
dt 



1 d 
2-Kr dr 



(8) 



Now 



we describe separately all the important ingredients entering this system of equations and 



all its necessary simplifications. 



2.1. Angular momentum flux. 

Goldreich & Tremaine (GT80) have considered the gravitational interaction of a gas disk with a 
planet embedded in it. They have demonstrated that a perturber on a circular orbit exerts a torque 
on the disk only in the immediate vicinity of the Lindblad resonances. Density waves launched at 
these locations carry angular momentum away from the satellite in the outer disk and towards 
it in the inner one. Lindblad resonances of order m ~ rp/hp ^ 1 located about one scaleheight 
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Table 1. Key to important symbols. 



Symbol 


Meaning 


Where defined 






o/^ n Q 4" 1 /~\Ti 1 /I 1 
cmJ.dliOll \^ J 




VloCUUo Ct-pdil cllliclcl 


89 ^ 


h 


QloK VciT/lCdji DCctlcIlclgllTj 


81 


7 


puiyLiupK.. lliClC X 


gz. 1 


T? 


total angular nioiiiciituin flux 


equation (9) 


II (n\ 

H'max ) 


luiiCTjioii Ciicxi acTjci TjOi C|Uc cuxon 


BCj^UailOIl \\y J 


-fj{r) 


QioijcLriCc-QcpciiQcriij ctnguicir iiioiriciixuiii nux 


cC[Ua.1jlOIl i^-lvJJ 


i 


QloljctllCc-ilKt? VctrictUic lOr Xilc bllOCK prop dig, ctlj lull 


cC[Ua.1jlOIl y^-^) 


(p[t) 


clllg,U.iai lIlUilKJilL Ulll LLclIlipiIl^ lUIlCLlOll 


t^qudirioii \^^) 


I 

Ip 


\ /I O fl 1 1 fYl 


equdLion irz) 


IVI I 


CiiLlCcll ^Idlit/Ldjiy lildjhs iUi d oLlUli^ WflVC llOillllitJclilLy 


equdLioii ^^o j 


1VI2 


ClillCcll jJlcLlicldiy lllcioo iUi ci bLlUli^ llilUl djLlUll iCCClUdCK 


equdiioii \i- j 


IVlf 


liOUCidii llidbb 


omip+inna ('\ \ kr 


X 


IdOldl LlibldllCt; ilUlil pidlltJb bCdicLl Uy hp 


equdiioli ^^J-Uy 


^sh 


dimensionless shocking distdnce 


equation (21) 


C 


proportionality constant in definition of Xsh 


equation (21) 


z 


radial distance from planet scaled by the shocking distance 


equation (22) 


a 


dimensionless surface density 


equation (19) 


PS 


dimensionless surface density at planet's location 


equation (19) 




drift velocity 


§2.2 & equation (26) 




relative difference in torques leading to migration 


equations (26) & (B4) 


V 


correction factor for drift velocity 


equation (26) 


Xs 


parameter characterizing strength of the feedback 


equation (27) 


zo 


cutoff distance 


equation (27) 


to 


characteristic timescale 


equation (29) 


K 


viscous parameter 


equation (30) 


At 


tidal parameter 


equation (31) 



Symbols with subscript "p" imply the value of the corresponding quantity at the planet's 
location. 
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~ hp from the planet are the strongest contributors to the angular momentum flux: lower order 

resonances are far from the planet and do not feel such strong tidal forcing while higher order ones 
are saturated because of the so-called "torque cutoff" (GT80). The total angular momentum flux 
across a cylinder of radius rp carried by all density wave harmonics is shown to be given by the 
expression (GT80) 

Fo = {GMpf^^ {^/^L.(Q)[2i^o(2/3) + Ki(2/3)]^} , (9) 

where denotes the modified Bessel function of order v and the function UmaxiQ) describes 
the strength of the torque cutoff. It depends only on the disk stability parameter Q mentioned 
previously in equation (4), and /J-max ~ 0.69 for disks with Q 3> 1. 

Beyond several hp from the planet, tidal perturbations behave basically like sound waves (Lin 
&; Shu 1968). Their nonlinear evolution finally leads to shock formation (GROl). After this happens 
the angular momentum flux of the density waves decreases because it gets transferred to the disk 
fluid. Taking this into effect account Rafikov (2001) has demonstrated that the angular momentum 
flux carried by the tidal perturbations Fj is given by 

23/2^3^ 

$(M„0 = / x\Mp,t,r,)dr, = (^^^ ' ^^^^ 



"PI/ 

Here A = (r/2)dJ7/dr is Oort's A constant, Ip is the Mach-1 distance [the distance from the planet 

at which shear velocity is equal to sound speed, Ip = (2/3) /ip for Kcplcrian disks], 7 is a polytropic 
index, and AIi is given by equation (3). The variable t{r) plays the role of the distance travelled by 
the wave from the planet, and the dimensionlcss flux $(Mi, t) was calculated in GROl for Mp = Mi. 
In equations (10)-(12) the surface density S and sound speed c can vary arbitrarily with radius r. 

It was also demonstrated in GROl that the density wave shocks when t reaches 

Ml 

tsh = 0.79-^. (13) 

Mp ^ ^ 

In this expression we have neglected an additional term equal to 1.89, which takes into account 
the finite distance required for the wake formation in the linear regime. One can do this only 
if Mp <^ Ml (which is always true in this paper), so that the region of linear wake formation is 
separated from the region of its weakly nonlinear evolution. Otherwise these regions overlap and 
the wave shocks before it is fully formed, significantly complicating the analysis. 

For t <tsh, the total angular momentum fiux of the waves Fj is conserved and has to coincide 
with Fo given by the expression (9); the value of $(t) for t < tgh is such that this condition is 



-8- 



fulfilled. Note also that although the polytropic index 7 enters equation (10) the dimensionless flux 
Fj does not actually depend on it because $ oc (7+1)^ (see GROl). Thus, we can rewrite equation 
(10) as 



The function ip{t) describes the damping after the shock is formed and is such that ip{t) = 1 when 
t ^ tgh, dip{t)/dt = when t = tgh, and (p{oo) = 0. Its behavior is shown in Fig. 7. In this paper 
instead of the true angular momentum flux dependence calculated by GROl we will be using for 
numerical convenience a simple analytical fit satisfying all of the above mentioned conditions and 
described in Appendix A. 

From the conservation of the total angular momentum (assuming a steady state) it follows 
that the corresponding angular momentum flux of the disk fluid is 



This expression together with equations (10)-(12) is to be substituted into equation (8) and then 
the surface density evolution could be calculated self-consistently because the dependence of F on 
S is taken explicitly into consideration by equation (12). 



Using equations (10)-(12) one can in principle study the temporal evolution of S taking into 
account the fact that Sq, cq, ^0, etc. vary on scales ~ r^. We prefer to study the problem 
in a simplified setting assuming that although gap formation is a nonlocal process it still occurs 
sufficiently close to the planet. Later we will determine the necessary conditions for this assumption 
to be true. 

This simplification allows us to set c = Cp, r = Vp, and to use the shearing sheet approximation 
to represent the background fluid flow. The only disk quantity which is assumed to vary in our 
consideration is the surface density S. By setting c = Cp wc imply that rapid variations of S do not 
violate the thermal balance in the disk that determines c(r). Here the disk temperature distribution 
is supposed to be the result of external irradiation by the central star and our assumption of constant 
c rests upon the absence of geometrical effects such as the shadowing by the gap edge (Lecar & 
Sasselov 1999; Dullemond et al. 2001). 

In this "quasilocal" approximation, introducing a new variable 



Fj{r) = Fo if , where 

ip{t) = ^Mp,t)/^Mp,tsh), 




(14) 




(15) 



2.2. Local approximation. 



(16) 
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we obtain from equation (12) that 

X 

t{x) = sign(a;)2-^/4 j 



.S(x') 



1/2 



\x'f''dx'. 



(17) 



If S(x) = Ep we reproduce the expression for t{x) obtained in GROl. 



The disk-planet interaction leads to the radial migration of the planet (GT80; Ward 1997a) 
with some drift velocity Vd- We consider the disk evolution in the coordinate system comoving with 
the planet. In this case 



S(r, t) —>Tj{r — Vdt, t) and 



dt 



dt ^'^ dr 



Usually it is found that migration is inward, i.e. Vd < (Ward 1986). Using the definitions 



a{r,t) = 



S(r,t) 
E(oo) ' 



PT. = 



E(oo)' 



(18) 



(19) 



one can rewrite equation (8) in the following form: 



+ 



92 



da Vd da _ 3u d^a 
dt Ip dx Ip dx"^ ' 7rE((X))Or2|p dx"^ 



(p sign(x) 



p^Mp }\xf^ 
2V4Miy (7V2 



dx' 



(20) 



Here E(oo) is the surface density of the disk far away from the planet; it is different from Ep in 
the case of viscous disks (r/ios 7^ 1, see §3.2). Then the ratio is not known a priori and has to 
be found from a self-consistent solution of the problem. We take E(oo) to be the same both inside 
and outside of the planet which means that global surface density gradients are neglected here as 
we have mentioned earlier. 

GROl has demonstrated that in the constant surface density disk, the wave shocks after prop- 
agating a distance (in x) 



Xsh ~ 1-4 



( l + l Mp 
V 12/5 Ml 



-2/5 



-2/5 



(21) 



away from the planet. In what follows we will always assume for simplicity that Q = 1.4 thus taking 
7 = 7/5 [see however Fridman &; Gor'kavyi (1999)]. In a disk with varying surface density Xsh is no 
longer the distance after which the shock forms. Still, it is convenient to change x to a new spatial 
coordinate z given by 



z 



r /M„ 



Ml 



2/5 



Xsh C 

In terms of z, the surface density in our particular problem varies on scales z 
condition for applying our "quasilocal" approach (see §2.2) is 



(22) 
Thus, the 



(23) 
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We check later in §4 if this is fulfihed in reahstic protoplanetary systems. 
With the aid of (22) equation (20) transforms to 

da__vd_ (MpV'^ da_3u_ f M^V^^ 
dt Qp \Mi) dz~ Clj \Mi) dz^ 



^1/2^5/2 } I /|3/2 

I{z) = sign(z)^^^^ J ^-^dz', Mf = ^{oo)h'. (25) 



Equation (24) does not have Mp inside the argument of (p and this permits one to extract an 
exphcit dependence on the planetary mass. 



2.3. Evolution equations. 

Equation (24) can only describe the surface density evolution when it is supplied with the 
information about the behavior of the migration velocity as a function of a. In Appendix B we 
demonstrate that 

"'^ = -iV^^"' "^^^^ - = 1-A.| / -/ l^^^-^rf., (26) 

is a correction factor for the drift velocity caused by the feedback from the surface density variations 
to the planetary migration; also 

3 r f M 

= -Wpr^ 7^ TT Zq = {UmaxXshY^- (27) 




PCtJ-maxhp \MiJ 

In the expression for v the migration feedback is represented by second term — it takes into account 
modification of the contribution from low-order Lindblad resonances, which is important because 
the planetary torque in the outer and inner parts of the disk due to high-order Lindblad resonances 
almost completely cancel. This means that even a small contribution from low-order resonances 
could produce a significant effect on the migration speed. When the feedback is absent the drift 
speed is determined only by global gradients of Sq and cq represented by factor /3 in equation (26) 
(see Appendix B). 

We will see later that in interesting cases the planet produces variations of a ~ 1 and a > 1 
in the direction of migration and cr < 1 in the opposite direction: a surface density increase forms 
in the disk in front of the planet, and a depression appears behind it. Thus density perturbations 
on both sides of the planet work in one direction — to enhance the effect of the disk in front of 
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the drifting planet and decrease it behind the planet [neglecting for a moment the fact that factor 
Py. in equation (26) could be different from unity]. They do not cancel each other (like in the case 
of migration due to variations of Sq) and produce 0(1) effect. This means that planetary-driven 
surface density perturbations tend to stall migration. However, these variations are strongest at 
distances ~ Xghhp S> hp from the planet, where torque generation is weak, thus they need to be 
large to compete with usual drift due to weak global gradients (for which cancellation effects are 
important) . 

Using this and the results of previous sections we can rewrite equations (24)- (25) to find the 
following form of the evolution equation : 

Pe^^ + = PE^A.g + A,^ (/(.))] , (28) 
r^l, to - X ^ lUlaAQmom) + i^i(2/3)]4"' ( ^V^' (29) 



V - « ^ {iA„(Q)[2if.(2/3) + (^1) ^i. (30) 

3 /MpV/^Mi 



^'^Acpywj M}- ^^^^ 

Here we have parameterized viscosity using the usual a-prescription: v = ahc and the correction 
factor for drift velocity v is given by equation (26). 

We will often use the numerical form of the parameters entering this system, namely 

A..0.16<j(*)'^ A..0.48^(*)"'\ (32) 

where we have used equation (4) and our calculation assumes that Q ^ 1, (3 ^ 7, and ^ ~ 1-4. 

Equation (28) describes the evolution of the disk surface density and its solution depends on 
5 dimensionless parameters: Xi,,Xt,Xs, to, and zq which themselves are combinations of hp/r-p, 
Mp/Mi, a, and Q. The parameter describes the effect of the planetary tidal torques on the 
surface density distribution, while tells us how strong the viscosity is. Obviously, for the gap to 
form the second term in the r.h.s. of (28) must be larger than the viscous one. The feedback which 
planetary migration receives from the surface density perturbations is described by the parameter 
As defined in equation (27) and it is strongest for large A^. Feedback also depends on the value of 
Zq [equation (27)] which represents the effect of the torque cutoff. 

To understand better the meaning of these parameters let us consider the time Ata which it 
takes for a planet to drift across a gap with the width Ar. Obviously, Ata = Ar/vd, where again 
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Vd is a drift velocity. Using equations (9) and (26) one can easily obtain that 

PhpMp ArQrpMi 

'"d^c—^-f, and Atd ^t^ti-- (34) 

Q Tp Ml c p hp Mp 

Time needed for viscous diffusion to fill the gap is clearly 

V Hp a \ hp J 

Finally, for the tidal torques to clear out a gap with width Ar an angular momentum H = 
Sil (rAr)^ must be supplied to the gas (see GT80). If we assume that angular momentum flux 
carried by the density waves is damped on the characteristic length ~ Ar, then the timescale for 
gap opening is 

In our particular problem the characteristic gap width and the wave damping distance are the 
same and given by Ar ~ hpXgh = (hp {Mp/Mi)~'^^^ [see equation (21)]. Then, substituting this into 
equation (34) we obtain that tg ~ ^td^ i-e- time to is a typical timescale for a gas drift through the 
gap caused by the planetary migration (in planet's reference frame). One can see that to ^ 
for Mp <C Ml which justifies our use of azimuthally averaged quantities if At ~ 1 [because to/\t is a 
characteristic timescale for the gap opening, see equation (28)]. Also for this choice of Ar one can 
easily find that At ~ At^/ Atopen, that is drift through the gap cannot replenish material repelled 
by planetary torques when At > 1, and gap could be opened if viscosity is absent. We will further 
confirm this conclusion in §3.1. Likewise, one can see that ^ Atd/ Aty which means that viscous 
diffusion is more important for the surface density evolution than effects of the migration when 
A;^ > 1. It is also clear that for tidal torques to overcome viscosity and open a gap in a viscous disk 
one needs Atopen ^ ^tv, or At > (see §3.2). 



3. Kinematic wave solution 

It is possible that for some values of parameters A,/, At, Xg solution for the surface density does 
not depend on time-like variable r. In this case, the profile of a is stationary in the reference 
frame migrating with the planet. Although planetary gravitational torques tend to repel the disk 
fluid, planet is mobile enough to migrate through the forming gap and stop its development. When 
such a time-independent solution is no longer possible, planet migration stops and surface density 
evolves with time to form a gap in the disk. 

To study such kinematic wave solutions we take da/dT = in equation (28). Integrating 
resulting equation with respect to z and using the fact that at infinity u = 1, da/dz = and 
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planetary torques are absent we obtain that 

v{a - 1) = p^'X.^ + A*^ b (/(z))] , I{z) = sign(z)^^^ j ^-^dz', (37) 



with V given by (26). 

WH89 used the same approach when studying gap formation with instantaneous density wave 
damping. In their case (p varies on scales ^ hp (wake generation requires a couple of disk scale 
lenghts to complete), i.e. dcp/dz ~ Xgh in our notation. Also, to significantly affect the migration 
speed one only needs surface density variations Aa hp/rp to be produced. Combined with the 
definition (31) this gives ~ {hp/rp)Mf as a value of their critical mass, in agreement with equation 
(1). In our case this mass is larger because wave damping occurs at greater distance from the 
planet. 

We now explore separately cases of purely inviscid disks and disks with nonzero viscosity. 



3.1. Discs with v = 

If the viscosity in the disk is absent equation (37) simplifies to 



a'^'ia - 1) = sign(z)^|A_A, \z\'/'^' {I{z)) , (38) 

where ip' (/) = dip (7) /dl. 

It follows from the properties of (p'{t) that cr = 1 for t < tgh, i-e. for \z\ < 1, meaning that 
Pe = 1 (see Appendix A). This implies that in inviscid disk nothing happens with the disk fluid 
until tidal wave shocks. After that, for \z\ > 1, p' < suggesting that for inward migration a > I 
if z < (inner disk) and o" < 1 if z > (outer disk) just like we mentioned before. The effect of 
this is to slow down the migration. We can rewrite (26) as 




(39) 



Since integrals are 0(1), one needs As > 1 for the drift velocity to be significantly affected. The 
solution of the system (38)- (39) depends on 2 parameters only, A^ and A^. 

In Fig. 1 we plot the surface density profiles in the vicinity of the planet for several values of 
Xf and Xs- Surface density increases in front of the migrating planet because there are two fiuxes 
of the disk fluid there — one caused by the migration, another by the planetary repulsion — and 
they are converging causing a surface density pileup. Behind the moving planet these fluxes are in 
the same direction — away from the planet. For large enough planetary mass repulsion removes 
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Fig. 1. — Surface density profiles for several values of At and Ag characterizing the strength of 
the planet's gravity and migration feedback. Surface density normalized by its value at infinity 
is plotted as a function of the distance from the perturber (normalized by the shocking distance 
Xshlp)- Black dot denotes the position of the planet; central star and direction of migration are to 
the left. 
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matter from the the disk behind the planet so rapidly that a corresponding surface density decrease 

cannot be replenished by the flux due to the planetary migration (this flux grows with Mp slower 
than the one due to the repulsion). As a result a gap is carved out in the disk behind the moving 
planet. 

To describe this process quantitatively let us notice that r.h.s. of equation (38) depends on a 
only through an integral over z. For 2; > r.h.s. of (38) varies from 

fmin = ^ min Ay^^ ^min (^^) 

to 0, where Zmin is the value of z for which ip' {I{z)) reaches minimum, and V'mm ~ —0.273 as 
described in Appendix A. At the same time l.h.s. of (38) varies from —2 x 3~^/^ to (it is 
minimized at t = 1/3). Thus, the solution of (38) does not exist if \fmin\ > 2 x 3~^/^ or if 

A > J_J^^!^ ,-3/2 .4^^ 
* q3/2a5/2|,^/ I "^min ■ V*"-) 
^ I rmin I 

Closer inspection of equation (38) reveals that a develops an infinite derivative when this happens. 
This would lead to the violation of the Rayleigh's criterion and instability associated with this will 
presumably help to clear out a gap on the outer side of the disk. 

Once the gap is formed, density waves launched in the outer disk arc not able to propagate to 
infinity and to transfer all their angular momentum to the disk material. Instead they refiect from 
the edge of a forming gap and return to the planet. Here they could interact gravitationally with 
planet thus canceling part (or all) of the wave angular momentum launched in the outer direction. 
It would diminish the influence of the outer disk on the migration and planetary drift will stall. 
For this to occur planet needs to absorb only a small fraction of one-sided angular momentum Fq, 
about hp/rp, which is obvious from the discussion in §2.3 and Appendix B. 

For As < 1 it turns out that Zmin ~ 1-2 and v ^ 1 (planet is still rapidly migrating, feedback 
effects are not strong enough to slow it down) when steady-state solution for a is no longer possible. 
This means that gap opens when 

At > 0.5 if Xs<l. (42) 

This corresponds to planetary mass 

Mp> Mt^ 2.3 Q-^/^ Ml if A, <1. (43) 

For stable disks {Q ^ 1) this criterion gives quite low mass at which gap forms although larger by 
~ {rp/hp)Q'^/'^ than the one suggested by Ward & Hourigan (WII89). 

For As > 1 the situation is somewhat different. Here the feedback from the density perturba- 
tions to the drift velocity of the planet becomes important. To study this region of the parameter 
space let us divide equation (38) by z^a^^"^ and integrate the result from —00 to —1 and from 1 to 
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Fig. 2. — Drift velocity of the planet of mass Mp with feedback Vdr, normalized by the drift velocity 
of the body with mass Ms [see equation (47)] without feedback, as a function of Mp/Ms- Velocity Vdr 
is shown by the thick solid line. Also shown are the drift velocity without feedback and analytical 
approximation for given by equations (44) &: (45) (both physical and unphysical roots are 
displayed). 
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oo, thus forming integrals present in r.h.s. of equation (39). Substituting them into (39) we find 
that 

v{l -v) = XsXtQiXs, Xt), where ^(A„ A*) = + 1 j ^T^j^^^- (44) 

It turns out that function £>(As, At) depends on A^, At only very weakly [which is reasonable because 
r.h.s. of (44) does not depend on the rapidly varying factor a — 1]. We confirmed this numerically 
and found that 

e(A„At) «0.31. (45) 

Since l.h.s. of equation (44) attains its maximum equal to 1/4 at v = 1/2, no solution of equation 
(44) for V exists when 

AsAt> , ,^ , . if A,>1. (46) 
MAs,At) 

Prom the condition (46) and definitions (27) & (31) we can find that the feedback stops 
migration when 

/ h \^/^^ 

Mp>M,^5.8 (^Q"'/J Ml if A, >1. (47) 

In Fig. 2 we show the dependence of the drift velocity on the mass of a perturber in the regime 
of strong migration feedback. Drift velocity is normalized by the drift speed which a body of mass 
Ms [equation (47)] would have if the feedback were absent (straight dashed line on this Figure). 
Analytical solution of equation (44) is also displayed and is in good agreement with numerical 
result. One can easily see that dvdr/dMp = — oo when Mp = Mg. This situation is analogous to 
what WH89 have found in their inertial limit calculations in the case of instantaneous damping. 
When a solution for the drift velocity does not exist any more planetary torques cannot support 
a constant flux of disk material seen in planet's frame and time dependent evolution commences 
leading to a gap formation and stalling the migration. Lin & Papaloizou (1986) and WH89 followed 
this process with time-dependent one-dimensional numerical calculations. 

The behavior of the drift velocity correction factor v is shown in Fig. 3. One can see the cutoff 
of the possible kinematic wave solutions caused by the migration feedback [predicted by equation 
(46)] on top of the plot where v rapidly diminishes. For small A^ the criterion given by (42) becomes 
important. The agreement between our simple analytical considerations and numerical results is 
very good. 

We can unite the gap formation conditions (43) and (47) in a single criterion on the limiting 
mass of the perturber: 



Mp . 

— - > mm 



Ml 



2.3 5.8 (q-'^) 



(48) 



Because Q ^ 1 we can immediately see from equation (48) that gap could be formed by a planet 
with a mass significantly smaller than Mi. 
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Fig. 3. — Contour plot of the correction factor v for the drift velocity of the planet [see equation (39)] 
, as a function of parameters At and As. Shaded region corresponds to the part of the parameter 
space where kinematic wave solution is not possible, migration stalls and gap opens. Vertical 
dotted line shows the restriction given by equation (42). Dashed line displays the limitation given 
by equation (46), important when the feedback is strong. One can see a very good agreement 
between the simple prescription given by equation (48) and actual numerical calculation. 
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Since our analysis has only exploited the most basic features of the damping function (p(t) which 
are well reproduced by our simple fitting expression (Al) we conclude that the results obtained 
here would be basically the same if one were to use the exact damping function calculated in GROl. 



3.2. Viscous discs. 

In the case of nonzero background viscosity present in the disk one faces a more complicated 
situation than in inviscid disk. Now the solution of kinematic wave equation (37) depends on 4 
parameters: At, A^, which we had before, X^,, which determines the effect of the viscosity on the 
solution, and Xgh, which affects drift velocity feedback. Parameter Xgh did not appear in the inviscid 
case because a was exactly 1 near the planet in that case. Viscosity changes this picture because 
now a varies near the planet, and this is important for the feedback, because tidal forcing is very 
strong in the immediate vicinity of the planet [see equation (B6)] . 

In the case u ^ one can integrate equation (37) once to obtain 



cr - 1 = PS 



At I vp^ 



cxp 



Xu 



{z - z') 



<p{l{z'))dz'-v{I{z)) 



(49) 



or 



(7-1 



PS 



A. 



ip ( I ( z + s 



-^{I{z)) 



exp(— s)(is. 



(50) 



If one integrates r.h.s. of (50) by parts and takes the limit Xi, — > 0, equation (38) is easily recovered. 

As we have mentioned before, factor in equation (50) is not equal to unity any more. It 
should be determined from the self-consistent solution of equation (50) instead. To solve this equa- 
tion we employ an iterative technique: at each step solution obtained from the previous iteration 
is substituted into the r.h.s. of (50) to obtain a better approximation until this process converges. 
As a first trial we use (j{z) = 1 and ps = 1- 

In Fig. 4a we display several surface density profiles for different X^ and fixed At = 0.52 (gap- 
forming value of At in the inviscid disk without the migration feedback) and A^ = (thus assuming 
that migration feedback is absent). For small Xi, profile is almost the same as in the inviscid disk; 
one can show that (cr — 1) ~ exp(— A~^) for small z. As Xi, grows and becomes comparable with At 
perturbations of a near 2 = become significant which has profound implications for the migration 
feedback if A^ ^ 0. However, as Xy > At viscous diffusion becomes so strong that surface density 
inhomogeneities in the vicinity of the planet become very small (as in the case Ajy = 3 in Fig. 4a). 

It turns out that in the viscous case (exactly like in the inviscid one) strong tidal perturbations 
can open a gap in the disk even without the drift velocity feedback. The presence of the viscosity 
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Fig. 4. — (top) Plot of the surface density as a function of distance from the perturber for different 

values of viscous parameter A,^: 0.01, 0.3, 1, 3. One can see how the viscosity smooths out density 
inhomogeneities in the disk and produces variations of a in the region \z\ < 1, leading to a stronger 
drift velocity feedback, (bottom) Maximum possible Aj for which planet is able to migrate without 
opening a gap as a function of A;^. It is assumed that feedback is absent: As = 0. One can see that 
viscosity tends to inhibit gap formation. 
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only changes the threshold value of the parameter so that it increases when X,y (and, correspond- 
ingly, viscosity) grows. The dependence of this limiting Aj on Ajy for Ag = (no feedback) is shown 
on Fig. 4b and is easy to understand: to open a gap the planet has to overcome not only the drift 
of the fresh disk material into a forming gap behind the planet but also the diffusion due to the 
viscous stresses which tends to smooth any inhomogeneities in a and fill the gap. Thus, larger Ajy 
require larger Xt to open a gap. 

When Xs > drift velocity feedback from the surface density perturbations could become 
important before At reaches the threshold value mentioned above. In Fig. 5 we show the boundary 
of the region in the At — Xg parameter space where a planet could migrate without opening a gap 
for several values of X^. As we have said before, for small Xg region where the gap cannot be 
cleared expands as Xi, increases confirming results for A^ = displayed in Fig. 4b. For larger Xg, 
however, the situation is more complicated: for a fixed strength of planetary torques (characterized 
by At) this region first rapidly contracts but then expands again with increasing X^- It could be 
explained as follows: when Ajy ^ At viscosity cannot strongly change the surface density near the 
planet, ps ~ 1 and drift velocity feedback is only slightly stronger than in the inviscid case. When 
viscosity grows so that X,^ ^ Xt, surface density perturbations are strong and viscosity is significant 
to sufficiently modify surface density in the planetary vicinity leading to a strong feedback: as 
it follows from equations (B6) and (26) even small deviations of a from a{0) = ps are strongly 
amplified by a factor ^ z"^ (and z ~ l/xgh ^ 1) in the expression for the feedback correction. 
This effect rapidly lowers the upper boundary of the region where planet could migrate and not 
open a gap. However, when Xi^ S> At strong viscosity smooths out all the inhomogeneities in a 
which planet tends to produce and this reduces the role of the feedback — for a fixed A^ critical A^ 
increases. Thus, growing viscosity tends to amplify migration feedback but then attenuates it as 
Xu exceeds At. 

In the viscous disk the absence of the steady-state solution of equation (50) caused by the 
feedback is not enough to open a gap, a fact emphasized in Lin &: Papaloizou (1986) and WH89. 
Planetary torques must also exceed viscous diffusion for this to happen. This additional requirement 
could be roughly described by a condition that 



where we have used equations (32)-(33). This is equivalent to saying that the time to open a gap 
by tidal torques must be shorter than the viscous diffusion time (sec §2.3). The condition (51) 
coincides with the corresponding criterion used by Lin & Papaloizou (1986) and WH89, but is 
different from the condition found in GT80 because only gap formation by a torque exerted at a 
single Lindblad resonance was considered in GT80. In Fig. 6 we plot the boundary of allowed 
region in the space of physical parameters: Mp/Mi and rp/hp for different values of dimensionless 
viscosity a and Toomre stability parameter Q. In depicting the boundary curves we have taken 
into account the condition (51) (straight line parts of the curves are due to this restriction). 

For a fixed viscosity, the planet can open a gap only if it is located to the right of the cor- 




(51) 
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Fig. 5. — Boundary of the region in which planet is able to migrate without opening a gap for 
different values of the viscous parameter \y in the Aj — plane. Planet migrates when At and A^ 
are small and gap is opened for large At and A^. It is assumed that Xgh = 5 is constant for all the 
planets (meaning that Mp/Mi is fixed). 
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responding curve in Fig. 6. For small Vp/hp (thick disks, small migration feedback) critical Mp 
is always larger than in the inviscid case and grows with increases both 

\u and Xi, grow according to equations (27) &: (30). This rapidly increases the role of the migra- 
tion feedback so that corresponding limiting Mp could become smaller than in the inviscid case. 
Complicated shape of the boundary curves for different values of a is due to the before mentioned 
complexity of the influence of the disk viscosity on the migration feedback. 

One can see from the Fig. 6 that for values of Vp/hp typical in protoplanetary disks {vp/hp ~ 
20—30) Mp <^ Ml if a is small (typically when a < 10~^). One can also notice that it is a migration 
feedback which initiates the evolution of the surface density leading to the gap formation once the 
condition (51) is fulfilled. 

Our consideration of the viscous case assumes that the only surface density gradients producing 
viscous flux are due to the planet-driven perturbations. However, background surface density 
varying on scales ~ rp also contributes to the viscous fluid flux. We neglected this effect completely 
when we employed our "quasilocal approximation" and omitted factors E/r leaving only dS / dr in 
our equations (see §2.2). If we were to include these factors in our analysis a picture similar to 
that described in Ward (1997a) would be found: for large viscosity there would not be an infinite 
derivative of the drift velocity with respect to Mp at some point as in the case shown in Fig. 2. 
Instead v^r would smoothly (but rapidly) decrease at some point as Mp grows and the migration 
would change gradually from type I to type II (Ward 1997a). We are unable to capture this effect 
working in the framework of our "quasilocal approximation" but this limitation cannot affect our 
major conclusions concerning the gap forming planetary mass. Also, when the viscosity is small or 
when the condition (23) is satisfied we should not worry about this effect at all. Using equations 
(32)-(33) we find that these additional viscous fiuxes could be neglected if Vp/hp » (Mp/Mi)~^/^ 
which is usually the case in realistic situations. 



4. Discussion 

4.1. Limitations of the analysis. 

The validity of our consideration is restricted not only by the requirements (23) and Mp <C Mi 
but also by the applicability of our nonlinear wave damping prescription. Its validity was extensively 
discussed in GROl and Rafikov (2001) and wc will not repeat this discussion here. For small 
enough viscosity {a < 10"*^) we found (sec Fig. 6) that Mp/Mi ~ 0.1. This means that typically 
tsh ~ 10 and Vsh ~ {2 — 2>) x hp so that the condition Xsh ^ 1 is marginally fulfllled (note also 
that tsh 3> io = 1-89 which is the offset value of t characterizing the wake generation region, see 
GROl). Thus, a wave typically shocks after travelling only several hp from the planet; assumption 
Xsh ^ Tp/hp is reasonable if Vp/hp > (10 — 20) which is the case in MMSN model described by 
equation (2) (vp/hp ~ 25 at 1 AU). In disks with higher viscosity, the wave will shock even closer to 
the planet (see §3.2) meaning that initially wake is not tightly wound, but this can hardly change 
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Fig. 6. — Boundary of the region in which planet is able to migrate without opening a gap for 

different values of dimensionless viscosity a and Toomre stability parameter Q in the Mp/Mi — 
—Vp/hp plane. We plot boundaries for Q = 10, 30, 60, 100 (labelled on the panels), and a = (solid 
line), a = 10~^ (long- dashed line), a = 10^"^ (short-dashed line), a = 10~^ (dotted line). Gap in 
the disk is opened when the system is to the right of the corresponding boundary. Straight line 
portions of the boundaries represent restriction > necessary for opening a gap in the strong 
drift velocity feedback limit. Complicated shape of the boundary curves for a 7^ is due to the 
complex interaction of the viscosity (affecting the surface density profile near the planet) and the 
drift velocity feedback. 
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our gap formation criterion. 

4.2. Time-dependent evolution 

After the planet grows so massive that a steady-state surface density distribution is no longer 
possible, time-dependent evolution commences. We do not study this process here. It is clear 
however that as soon as the gap opening criterion is fulfilled and disk viscosity is overcome by 
planetary torques, gap clearing starts behind the planet (in the part of the disk opposite to the 
direction of migration) . As we have mentioned before this makes propagation of the density waves 
in the outer disk impossible, and planet migration which is the result of the small difference of 
the torques on both sides of the disk, has to stall. When this happens the planet gradually repels 
the surrounding material, carving out a gap in both parts of the disk. Far from the planet gap 
expansion can be stopped by the disk viscosity or tidal torques from other planets in the system. 

This general picture of gap formation (after the corresponding gap-opening criterion is fulfilled) 
was confirmed by the one-dimensional time-dependent numerical calculations in Lin & Papaloizou 
(1986) and WH89. Of course they used a different gap-forming criterion because of their assumed 
instantaneous damping of the planet-generated density waves. But because the reason for the 
gap formation in our case is basically the same as in theirs — migration feedback and strong tidal 
torques exceed the viscous spreading — we expect the time-dependent transition from type I to type 
II migration in our setting to be similar to what these authors have found in the local damping 
approximation. The only differences would be the critical mass at which this happens and the 
transition timescale. Our results for the gap-opening criterion also seem to be in general agreement 
with the time-dependent two-dimensional simulations carried out by LP93, Bryden et al. (1999), 
and Nelson et al. (2000), although it is hard to make a direct quantitative comparison because our 
criterion requires knowledge of a larger number of parameters than is quoted in these studies. 

Typical timescale for a gap to form is to/A*. Using equation (29) or (33), assuming Mp = O.lMi 
and At ~ 1, one can find that for MMSN parameters this time is about 5 x 10^ 7 x 10^ 

yr at 1 AU and ~ 2 x 10^ ^ 4 x 10^ yr at 5 AU (scmimajor axis of Jupiter's orbit). Such 
long timescales can put some restrictions on the numerical schemes which would be able to check 
the validity of our gap-opening criterion (they should also have high spatial resolution, typically a 
fraction of hp, and be able to follow the formation and evolution of weak shocks in a disk, although 
numerical viscosity would be a primary concern). 

4.3. Applications 

Using the theory developed in §3 we can determine the planetary mass at a particular location 
in a disk for which gap formation is expected and migration switches to a type II mode. For 
the MMSN parameters represented by equation (2) one obtains Vp/hp 25, Q ~ 70, and Mi ^ 
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14 Mq at 1 AU. Using Fig. 6 we find that for a < 10" a gap is opened by a planet with 
Mp ~ (0.12 — 0.15)Mi ~ (1.5 — 2) M^. Thus, if the Earth was immersed in a gaseous disk at 
the end of its formation, it probably could not open a gap (nonzero disk viscosity only strengthens 
this conclusion). Then type I migration would have caused it to drift to the Sun on a timescale 
< 10^ yr, but this apparently did not happen. This result supports the usual view according to 
which the final accumulation of terrestrial planets via the runaway coagulation of planetesimals 
or giant impacts of Moon-sized planetary embryos occurred after the nebula was dispersed in the 
inner Solar System. Then Earth would not have migrated at all and the question of its survival 
would not arise. 

At 5 AU one finds that rp/hp 16, Q sa 45, and Mi Ri 50 M®; then critical Mp/Mi 
(0.15 — 0.2)Mi ?s (7 — 9)M0 for a < 10~^. Thus, even the usual rocky core of Jupiter alone 
(without gaseous envelope) would likely be able to open a gap; this would significantly slow down 
its migration towards the Sun and leave Jupiter on its orbit far from the central star. For some 
reason this clearly did not happen in systems harboring extrasolar giant planets close to their 
parent stars, and this issue apparently deserves further study. The value of the critical mass also 
raises another question: how has Jupiter managed to acquire its huge gaseous mass? It might be 
that the timescale for the gap formation at the Jupiter's location is long enough (see §4.2) for the 
planet to accrete all its mass during the gap-opening stage if the core instability (Mizuno 1980) was 
operating since the very beginning (core instability and associated planetary mass growth might 
have actually triggered the gap formation). But obviously more detailed consideration of the gas 
accretion process is needed to definitely answer this question. 

An interesting result following from this analysis is that when a gap is opened in an inviscid 
disk, there is still some material remaining at the planet's orbit. This material has the form of a 
ribbon with radial width equal to 2hpXsh because waves launched by the planet cannot shock prior 
to travelling a minimum necessary distance, and, thus, cannot transfer their angular momentum to 
the disk and cause its evolution. The effects of wave action refiection from the edge of the ribbon 
or the presence of the moderate viscosity in the disk might lead to the dissipation of such a gas 
torus. But if these effects are not very strong then this ribbon phenomenon might have observable 
manifestations. 

5. Conclusions 

Using a realistic damping prescription for tidally-induccd density waves wc have studied the 
conditions necessary for a gap to be formed in a gas disk in the vicinity of a planet. It was shown 
that for small enough planetary mass a steady-state solution for the surface density perturbations 
exists in the reference frame migrating with the planet. Then the details of the surface density 
distribution in the disk depend on only 4 parameters: the aspect ratio in the disk hp/vp, the viscosity 
(represented here by its dimensionless analog a), the Toomre stability parameter of the disk Q, and 
the ratio of the planetary mass to a fiducial mass Mi defined in equation (3). Only systems with 
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a < 10"^ were considered because larger viscosity could violate our damping prescription based 

on the nonlinear wave dissipation. We have demonstrated that in disks with Q between ~ 10 and 
~ 100 (which should be typical in passive protoplanetary disks irradiated by their central stars) 
a gap is opened when the planetary mass reaches (1 — 10)M^, depending on the disk viscosity 
and the planet's location in the nebula. Planets further away from the central star must be more 
massive to repel the gas in their vicinity. We obtained an analytical criterion for a gap-forming 
planetary mass in inviscid disks [sec equation (48)]. Our result for this critical mass is in between 
the prcvioTis estimates of this quantity obtained by Ward & Hourigan (1989) and Lin & Papaloizou 
(1993) because of the different density wave damping function (WH89) and moderate requirements 
for the nonlinearity of the wave (LP93) . 

The apparatus developed here for following the disk surface density evolution and studying 
stationary structures in the disk such as the kinematic wave solutions (see §3) could be generalized 
to other astrophysical problems; for example it could be used to study global disk evolution or to 
extend our gap formation analysis to a system containing several planets. 
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A. Angular momentum flux if{t). 

GROl have calculated the dependence of the angular momentum flux of shock damped density 
waves using the full numerical solution for the propagation and damping of the wake. Here for 
simplicity we use a simple matching function which reproduces the main features of the numerical 
solution: 

0, t < tsh, 

The behavior of the true solution and the matching function is shown in Fig. 7. For a planet- 
generated wake profile tsh = 0.79 (neglecting the fact that there is some distance to over which 
the wake profile forms). Our simple fit never deviates from the true ip{t) by more than 17%. 
Asymptotically (p{t) oc t~^/'^ for t —>■ oo, as in the numerical solution. Also d(p{t)/dt = when 

t = tsh- 

To study the surface density profile near the planet in inviscid disks it is important to know 
the behavior of (p'{t) = dip{t)/dt. For our assumed matching function (Al) one can find that 

(p'{t) oc {tsh -t), t- tsh < 1, -t"^/^, t > 1, and (p'{t) = for i < tsh- (A2) 

The function (p'{t) reaches its minimum value equal to ^'^^^ = —0.273 at tmin = 1.43. 



t > tsh- 



B. Drift velocity 



From the global conservation of the angular momentum we find the velocity of migration driven 
by the planetary gravitational torques: 

dF 



MprpQp J 



dr 



dr, 



(Bl) 



where the integral is taken over the whole disk. One can write dF/dr = 'E{r)df /dr, where F is the 
angular momentum flux produced by the planet per unit surface density at distance r. GT80 have 



for 



fpl Z hp- 



demonstrated that function df/dr rapidly falls off like |r 

Ward (1986) has shown that df/dr = sign{r—rp)df /dr\o[l+0{hp/rp){r—rp)/hp], where df/dr\o 
is an even function of r — r^ and terms of order hp/rp arise because of the asymmetries in the torque 
generation pattern intrinsic for Keplerian disks. The surface density itself may vary on different 
scales. In our particular case we are interested in the variations on scales of order rp [represented 
by So(r)] which lead to migration in the first place, and variations on scales ~ hpXgh rp, which 
result from the planet-driven surface density evolution. In the vicinity of the planet S could be 
written in the following form: 

S-So(r)" 



S(r) = Eo(r) 



1 + 



So(r) 



^^2hpdlnEo^ 
3rp dlnr 







X 





(B2) 
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Fig. 7. — (top) Behavior of the dimensionless angular momentum flux (p (solid line) and a simple 
fit (dashed line) given by equation (Al). The flux diminishes as t increases due to shock damping. 
{bottom) Plot of ip'{t) calculated using the fitting function (Al). 
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{p-^(T — l)dx 



(B3) 



The replacement of T,o{r) with Ep in the last bracket assumes that variations due to planetary 
torques are more rapid than background ones (because damping length is supposed to be much 
shorter than Vp). 

Substituting this in (Bl) and recalling that df/dr\o falls off beyond several hp from the planet 
one finds that 

f dF , ^ \hp Sp /• , df 

Here /3 is a constant of the order of unity which depends on the global gradients of the surface 
density k = — dlnS/dlnr and temperature / = — dlnT/dlnr in the disk. The first term in brackets 
arises as a result of the asymmetry in torques produced in the inner and outer parts of the disk. 
Ward (1986) advocates that for Q = oo {Q \s a, Toomre stability parameter) 

/3 X 2 |^/x^„,(Q)[2Ko(2/3) + Ki(2/3)]2| = 6.5(1 + 0.06A; + 1.2/), (B4) 

and for Q = 2 r.h.s. of (B4) changes to 72(1 — 0.19A; + 0.95Z). It means that the surface density 
decreasing from the center of the disk leads to a very weak acceleration of drift in Q = oo disk and 
slows it down for moderate values of Q. Further we will assume for simplicity that the expression 
in parentheses in r.h.s. of (B4) is equal to 2 (which corresponds to /3 « 7). 

Using (B3) we obtain that 



WFq hp iTp^p [ df 

lVJ.pl p^Lp lp j^J ' U J LL^ 



{p^'^a - l)dx (B5) 



is a correction factor for the drift velocity caused by the feedback from the surface density variations 
to the planetary migration. 



For large x = {r — rp)/lp it could be demonstrated (GT80) that 

dx 



^ (B6) 



y ,,3 ^4 ■ 

r-max-^ 



For X < 1 torque rapidly decreases ("torque cutoff"). For simplicity we will assume that the 
dependence given by equation (B6) holds true even for x ~ 1, then we only need to correct the 
cutoff value of x so that the integral of (B6) gives us the right amount of the angular momentum 
flux Fq at infinity. It is done by setting df/dx\o = for \x\ < p^ax'i ^r \x\ > p^ax '^^ assume 
df /dx\o to be given by equation (B6). Using equation (B5) we finally obtain equation (26) and 
(27). 

Ideally, one should also take into account gradient of S on short damping scale in the expression 
for df /dr as it was done by Ward (1986). This could produce some modification because surface 
density gradients displace the positions of the Lindblad resonances and modify the torque cutoff. 
However, while variations of S produce contribution to the drift oc Act, presence of the gradients 
of S gives rise to the effects of order da/dx ~ Aa/xgh- This implies that contribution due to the 
gradient of the surface density is ~ 1/ Xgh compared with the last term in (B3) and is unimportant 
for Mp < Ml. 



